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1.  INTRODUCTION 


The  Boston  College  Institute  for  Scientific  Research  has  performed  innovative  research 
and  analysis  to  increase  the  knowledge  of  space  effects  on  Air  Force  systems.  Sections  in  this 
report  address  a  wide  range  of  efforts  from  applying  chemical  techniques  to  the  development 
of  new  technologies  for  Air  Force  space-based  systems  to  developing  a  database  of  various 
measurements  made  by  the  Defense  Meteorological  Satellite  Program  (DMSP).  In  addition, 
our  research  explored  solar  and  interplanetary  disturbances  and  their  coupling  to  geomagnetic 
storms.  Another  study  involved  research  and  development  of  highly  miniaturized  satellite 
skin  sensor  technology  capable  of  measuring  ionospheric  neutral  particles.  This  real-time,  in- 
situ  capability  enables  distributed,  remote  density  measurements  via  small  satellites  including 
the  Air  Force  C/NOFS  (Communications/Navigation  Outage  Forecast  System)  satellite 
mission. 

All  aspects  of  the  research  process  were  included  in  our  efforts.  From  data  acquisition, 
analysis,  theoretical  explanation,  database  development  to  summary  of  our  efforts.  Much  of 
our  work  was  published  in  peer  reviewed  journals  as  evidenced  in  Section  7  with  more  than 
45  publications  over  the  course  of  this  contract. 


2.  CHEMICAL  PROCESSES 

Research  performed  under  this  contract  included  experimental  and  theoretical  studies 
conducted  in  the  plasma  chemistry  and  mass  spectroscopy  laboratories,  located  on-site  at 
AFRL/RVBXT.  Experiments  included:  the  selected  ion  flow  drift  tube,  the  high  temperature 
flowing  afterglow,  the  turbulent  ion  flow  tube  and  the  flowing  afterglow  Langmuir  probe. 
These  are  fast  flow  kinetics  devices  for  measuring  ion  and  electron  kinetics  over  a  wide  rage 
of  temperatures  and  pressures. 

The  overall  purpose  of  this  research  was  to  understand  and  mitigate  the  fundamental 
processes  that  affect  a  wide-variety  of  Air  Force  interests  where  plasmas  play  a  prominent 
role,  e.g.  the  ionosphere,  reentry  vehicles  and  combustion. 

Data  measured  included  ion-molecule  reaction  rate  constants  and  branching  ratios  over 
extended  temperatures  and  pressure  ranges.  Research  efforts  also  included  engineering  and 
maintenance  tasks  to  design,  assemble  and  test  the  experimental  apparatuses.  This  required 
routine  cleaning,  maintenance  of  the  vacuum  pumps,  as  well  as,  leak  checking,  plumbing  gas 
lines,  and  repairing  electronics.  Instrument  upgrades  were  performed  as  necessary.  Near  the 
initial  phase  of  this  contract,  a  new  custom-made  heating  system  was  installed  on  the 
turbulent  ion  flow  tube  (TIFT)  to  ensure  uniform  heating  over  the  length  of  the  flow  tube 
including  better  temperature  measurement  and  monitoring. 

This  work  also  included  the  development,  validation  and  application  of  modeling  and 
simulation  tools  to  predict  plasma  formation  on  Air  Force  reentry  and  hypersonic  vehicles  and 
to  specify  the  impact  of  plasma  on  reentry  signatures  and  communication/navigations 
systems. 

All  aspects  of  the  research  were  supported  including  data  acquisition,  data  analysis, 
equipment  maintenance  and  upgrade  designs.  All  data  sets  were  delivered  to  AFRL  as  they 
became  available.  Results  of  our  research  were  presented  at  national  and  international 
conferences  and  published  in  peer  reviewed  journals.  Section  6.0  lists  nearly  40  publications 
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that  resulted  from  our  efforts.  Some  of  the  more  notable  research  results  are  summarized  in 
the  following  sections. 

2.1.  Kinetics  Studies  of  a  Variety  of  Positive  and  Negative  Ions 

Kinetics  studies  of  a  variety  of  positive  and  negative  ions  reacting  with  the  GX 
surrogate,  dimethyl  methylphosphonate  (DMMP),  were  performed.  All  protonated  species 
reacted  rapidly,  that  is,  at  the  collision  limit.  The  protonated  reactant  ions  created  from 
neutrals  with  proton  affinities  (PAs)  less  than  or  equal  to  the  PA  for  ammonia  reacted 
exclusively  by  nondissociative  proton  transfer.  Hydrated  HsO^  ions  also  reacted  rapidly  by 
proton  transfer,  with  25%  of  the  products  from  the  second  hydrate,  H30^(H20)2,  forming  the 
hydrated  form  of  protonated  DMMP.  Both  methylamine  and  triethylamine  reacted  exclusively 
by  clustering.  NO^  also  clustered  with  DMMP  at  about  70%  of  the  collision  rate  constant. 
and  ©2^  formed  a  variety  of  products  in  reactions  with  DMMP,  with  02^  forming  the 
nondissociative  charge  transfer  product  about  50%  of  the  time.  On  the  other  hand,  many 
negative  ions  were  less  reactive,  particularly,  SFs”,  SFe”,  COs”,  and  NOs”.  However,  F~,  O", 
and  02~  all  reacted  rapidly  to  generate  m/z  =  1 09  amu  anions  (P03C2H6~).  In  addition,  product 
ions  with  m/z  =  122  amu  from  H2^  loss  to  form  H2O  were  the  dominant  ions  produced  in  the 
0~  reaction.  N02~  underwent  a  slow  association  reaction  with  DMMP  at  0.4  Torr.  G3(MP2) 
calculations  of  the  ion  energetics  properties  of  DMMP,  sarin,  and  soman  were  also  performed. 
The  calculated  ionization  potentials,  proton  affinities,  and  fluoride  affinities  were  consistent 
with  the  trends  in  the  measured  kinetics  and  product  ion  branching  ratios.  The  experimental 
results  coupled  with  the  calculated  ion  energetics  helped  to  predict  which  ion  chemistry  would 
be  most  useful  for  trace  detection  of  the  actual  chemical  agents. 

2.2.  Kinetics  of  Ion-Molecule  Reactions 

The  rate  constants  and  product  ion  branching  ratios  have  been  measured  in  a  selected 
ion  flow  tube  (SIFT)  at  298  K  for  a  variety  of  positive  and  negative  ions  reacting  with  2- 
chloroethyl  ethyl  sulfide  (2-CEES),  a  surrogate  for  mustard  gas  (HD).  This  series  of 
experiments  is  designed  to  elucidate  ion-molecule  reactions  that  have  large  rate  constants  and 
produce  unique  product  ions  to  guide  the  development  of  chemical  ionization  mass 
spectrometry  (CIMS)  detection  methods  for  the  chemical  weapon  agent  using  the  surrogate 
instead.  The  negative  ions  typically  used  in  CIMS  instruments  are  essentially  unreactive  with 
2-CEES,  that  is,  SFe”,  SF4~,  CF30~,  and  C03~.  A  few  negative  ions  such  as  N02~  and  N03~ 
undergo  three-body  association  to  give  a  unique  product  ion,  but  the  bimolecular  rate 
constants  are  small  in  the  SIFT.  Positive  ions  typically  react  at  the  collisional  limit,  primarily 
by  charge  and  proton  transfer,  some  of  which  is  dissociative.  For  ions  with  high  proton 
binding  energies,  association  with  2-CEES  has  also  been  observed.  Many  of  these  reactions 
produced  ions  with  the  2-CEES  intact,  including  the  parent  cation,  the  protonated  cation,  and 
clusters.  G3(MP2)  calculations  of  the  thermochemical  properties  for  2-CEES  and  mustard 
have  been  performed,  along  with  calculations  of  the  structures  for  the  observed  product 
cations.  Reacting  a  series  of  protonated  neutral  molecules  with  2-CEES  brackets  the  proton 
affinity  (PA)  to  between  812  ((CH3)2CO)  and  854  (NH3)  kJ  mol"'.  G3(MP2)  calculations  give 
a  PA  for  2-CEES  of  823  kJ  moP'  and  a  PA  for  mustard  of  796  kJ  moP',  indicating  that  the 
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present  results  for  2-  CEES  should  be  directly  transferable  to  mustard  to  design  a  CIMS 
detection  scheme. 

2.3.  Reaction  Kinetics  of  P02Cr,  PO2CI2',  POCI2',  and  POCb*  with  O2  and  O3 

Rate  constants  and  product  ion  branching  fractions  for  the  gas-phase  reactions  of  O2  and 
O3  with  the  anions  (a)  P02Cr,  (b)  POCI3’,  (c)  POCI2',  and  (d)  PO2CI2'  were  measured  in  a 
selected-ion  flow  tube  (SIFT).  The  kinetics  were  measured  at  temperatures  of  1 63-400  K  and 
a  He  pressure  of  0.4  Torr.  Only  P02Cr  reacts  with  O2  to  a  measurable  extent,  having 
^163-400  K)  =  1.1  X  10'*(77K)'''®  cm^  molecule''  s'',  while  O3  reacts  with  all  of  the  anions 
except  PO2CI2’.  The  fitted  rate  constant  expressions  for  the  O3  reaction  with  anions  a-c  are  as 
follows:  /ta(163-400  K)  =  3.5  x  lO-^CTTK)''  ^  /tb(163-400  K)  =  4.0  x  lO'^Cr/K)'' ^  and 
^c(163-400  K)  =  3.7  x  10'’(77K)''  ''  cm^  molecule''  s''.  Calculations  were  performed  at  the  G3 
level  of  theory  to  obtain  optimized  geometries,  energies,  and  electron  affinities  (EAs)  of  the 
reactant  and  product  species,  as  well  as  to  determine  the  reaction  thermochemistry  to  help 
understand  the  experimental  results.  The  PO^Cl/  anions  that  have  lower  electron  binding 
energies  (eBE)  and  higher  spin  multiplicities  are  more  reactive.  The  doublets  are  more  labile 
than  the  singlets.  How  the  extra  electron  density  is  distributed  in  the  anion  does  not  predict  the 
observed  reactivity  of  the  ion.  The  reactions  of  P02Cr  with  O2  and  O3  yield  predominantly 
PO3'  and  PO4'.  The  reaction  of  POCI2’  with  O3  yields  mostly  Cl'  and  PO2CI2',  while  the 
POCI3'  reaction  with  O3  yields  mostly  O3'  and  PO2CI2’. 

2.4.  Electron  Attachment  to  SF5X  Compounds 

Rate  constants  and  ion  product  channels  have  been  measured  for  electron  attachment  to 
four  SF5  compounds,  SF5C6H5,  SF5C2H3,  S2F10,  and  SFsBr,  and  these  data  are  compared  to 
earlier  results  for  SFe,  SF5CI,  and  SF5CF3.  The  present  rate  constants  range  over  a  factor  of 
600  in  magnitude.  Rate  constants  measured  in  this  work  at  300  K  are  9.9  ±  3.0  x  10'* 
(SF5C6H5),  7.3  ±  1.8  X  10'^  (SF5C2H3),  6.5  ±  2.5  x  10''°  (S2F10),  and  3.8  ±  2.0  x  10''°  (SFjBr), 
all  in  cm^  s''  units.  SF5'  was  the  sole  ionic  product  observed  for  300-550  K,  though  in  the 
case  of  S2F10  it  cannot  be  ascertained  whether  the  minor  SF4'  and  SFb'  products  observed  in 
the  mass  spectra  are  due  to  attachment  to  S2F10  or  to  impurities.  G3(MP2)  electronic  structure 
calculations  (G2  for  SFsBr)  have  been  carried  out  for  the  neutrals  and  anions  of  these  species, 
primarily  to  determine  electron  affinities  and  the  energetics  of  possible  attachment  reaction 
channels.  Electron  affinities  were  calculated  to  be  0.88  (SF5C6H5),  0.70  (SF5C2H3),  2.95 
(S2F10),  and  2.73  eV  (SFsBr).  An  anticorrelation  is  found  for  the  Arrhenius  .4-factor  with 
exothermicity  for  SF5'  production  for  the  seven  molecules  listed  above.  The  Arrhenius 
activation  energy  was  found  to  be  anticorrelated  with  the  bond  strength  of  the  parent  ion. 

2.5.  Acidity  of  a  Nucleotide  Base 

Experiment  and  calculations  are  used  to  show  that  the  gas-phase  acidity  of  uracil  is 
comparable  to  that  of  HCl.  The  gas-phase  acidity  of  uracil  (denoted  here  by  U)  was  bracketed 
by  proton-transfer  measurements  involving  U  and  various  reference  acids  (denoted  here  by  A) 
of  known  gas-phase  acidity.  Rate  constants  for  proton  transfer  from  the  reference  acid  A  to 
the  conjugate  anion  of  uracil,  (U-H)',  were  measured  in  a  selected-ion  flow  tube  at  298  K. 
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Rate  constants  for  proton  transfer  from  U  to  ions  (A-H)'  were  measured  at  467  K  in  a  flowing- 
afterglow  Langmuir  probe  apparatus.  Here,  (U-H)  or  (A-H)  indicates  a  U  or  A  molecule 
which  is  missing  an  H  atom,  respectively.  The  result  is  A/7°acid(uracil)  =  333  ±  5  kcal  mol'' 
and  AG°acid(uracil)  =  326  ±  5  kcal  mol''  at  298  K,  which  agrees  with  earlier  work.  Thermal 
electron  attachment  to  uracil  was  found  to  be  too  slow  to  permit  measurement  of  a  rate 
constant,  consistent  with  the  gas-phase  acidity  given  above.  G3  and  G3(MP2)  calculations  are 
reported  for  uracil,  and  for  the  each  of  the  (U-H)  radicals  and  (U-H)'  ions  that  result  from  H  or 
H^  loss  from  each  of  the  four  hydrogen  sites  of  U  (on  the  Nl,  N3,  C5,  and  C6  positions). 

From  the  calculated  total  energies  we  obtain  the  gas-phase  acidity  of  uracil,  the  four  U-H 
homolytic  bond  strengths,  and  the  electron  affinities  of  the  four  possible  fragment  radicals. 

We  confirm  earlier  work  that  the  most  acidic  site  in  uracil  is  at  the  Nl  site;  this  site  is  where 
uracil  becomes  covalently  bonded  to  a  carbon  of  the  ribose  sugar  in  RNA.  G3  calculations  for 
the  Nl  site  at  298  K  give  A/7°acid(uracil)  =  334.5  kcal  mol''  and  AG°acid(uracil)  =  327.1  kcal 
mol''  at  298  K,  in  good  agreement  with  the  experiment.  The  weakest  H-atom  bond  enthalpy 
(at  the  Nl  site)  is  calculated  to  be  101.8  kcal  mol''. 

2.6.  Electron  Attachment  to  Ch  from  300  to  1100  K 

Rate  constants  for  dissociative  electron  attachment  to  CI2  have  been  measured  from  300 
to  1 1 00  K  in  a  high-temperature  flowing-afterglow  Langmuir-probe  apparatus.  R-matrix 
calculations  have  been  carried  out  which  compare  well  with  the  present  measurements  with 
possible  deviation  at  the  highest  temperatures.  The  attachment  rate  constants  do  not  show 
Arrhenius  behavior.  The  temperature  dependence  of  the  calculated  rate  constants  for 
successive  vibrational  levels  provides  insight  as  to  this  behavior.  While  the  lowest  vibrational 
level  of  CI2  dominates  attachment  at  low  temperatures,  the  rate  constant  is  not  flat  with 
temperature  because  of  the  /7-wavecharacter  of  the  attachment  process.  The  non-Arrhenius 
behavior  is  due  to  a  conflict  between  the  increase  in  attachment  cross  section  with  vibrational 
level  (temperature)  and  the  decline  in  the  crosss  section  with  electron  energy  above  50meV. 

2.7.  Electron  Attachment  and  Detachment 

Measurements  are  reported  of  rate  constants  for  electron  attachment  to  C6F5X 
(X=Cl,Br,l)  and  thermal  electron  detachment  from  C6F5C1-  over  the  temperature  range  300- 
550  K  in  133  Pa  of  He  gas  in  a  flowing-afterglow  Langmuir-probe  apparatus.  This  is  the  first 
case  we  know  of  where  the  parent  anion  has  sufficiently  low  electron  detachment  energy  that 
detachment  (from  C6F5C1-  in  this  case)  has  been  observed  in  competition  with  a  channel  for 
dissociative  electron  attachment  yielding  a  thermally  stable  anion  (here,  C1-).  Because  of  this 
competition,  it  is  shown  that  a  simple  mass  spectrometric  determination  of  the  product 
branching  fractions  at  long  times  will  lead  to  erroneous  results  at  elevated  temperatures.  The 
electron  density  profiles  provide  evidence  for  a  new  plasma  decay  process  involving  the 
detaching  and  nondetaching  anions  trapped  in  the  space  charge  field  of  the  positive  ions. 
Electron  attachment  rate  constants  were  found  to  be  1 .0^  1 0-7.  1 . 1  x  1 0-7.  and  2. Ox  1 0-7  cm3 
s-1.  at  300  K.,  for  C6F5C1,  C6F5Br,  and  C6F5I,  respectively,  estimated  accurate  to  ±25% 
except  for  C6F51,  where  there  is  ±30%  uncertainty.  Rate  constants  for  C6F5C1  changed  little 
over  our  temperature  range,  while  those  for  C6F5Br,  and  C6F5I  increased  with  temperature. 
Electron  detachment  occurred  only  for  C6F5C1-  in  our  temperature  range.  Detachment  rate 
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constants  were  immeasurable  at  room  temperature  but  approached  4000  s-l  at  550  K.  From 
these  data  the  electron  affinity  (EA)  for  C6F5C1  was  determined,  EA  (C6F5C1)=0.75±0.08 
eV.  G3(MP2)  calculations  (based  on  Moller-Plesset  perturbation  theory)  were  carried  out  for 
the  neutral  and  anion  and  yielded  EA(C6F5C1)=0.728  eV. 

2.8.  Strong  Emission  from  CaO  in  the  500  to  900  nm  Region 

Research  revealed  some  interesting  spectral  data  which  showed  strong  emission  from 
CaO  in  the  500  to  900  nm  region.  The  purpose  of  this  section  is  to  document  the  procedure 
that  was  used  to  make  a  rough  prediction  of  these  spectra. 

The  spectra  of  particular  interest  are  generated  by  the  vibration-rotation  transitions  from 
the  a'Z  electronic  state  to  the  X'l  ground  electronic  state.  Because  the  anharmonicity  of  the 
vibrational  sublevels  of  these  electronic  states  is  relatively  small,  these  spectra  are  grouped 
into  bands  corresponding  to  the  vibrational  quantum  number  change.'  The  anharmoinicity 
leads  to  a  “blurring”  of  the  band  heads.  The  rotational  spreading  of  the  bands  is  toward  the 
red. 

An  example  of  these  spectra  is  shown  for  a  re-entry  object  in  Figure  1 .  The  peaks  at  the 
right  of  the  figure  are  labeled  with  the  change  in  vibration  quantum  number  during  the 
transition.  Peaks  are  easily  seen  for  Av  =  0,  1,2  and  3.  There  are  weaker  transitions  from  a 
lower  lying.  A'  P,  electronic  state  to  the  ground  state  which  complicate  this  spectrum  ’  under 
some  conditions,  but  they  are  not  of  interest  to  this  study.  The  spectra  labeled  green  and 
orange  a  complex  spectra  originating  primarily  from  higher  electronic  states  of  CaO.'*'’ 

The  dashed  line  in  the  figure  is  representative  of  a  black  body  curve  at  1750  K.  This  is 
representative  of  a  solid  surface  radiating  at  the  temperature.  The  spectrum  has  been 
transmitted  through  a  large  atmospheric  path  and  thus  shows  the  oxygen  “atmospheric”  band 
absorption  at  7600  and  6900  A.  The  usual  sodium  radiation  at  5900  A  is  also  expected. 

A  computation  of  the  spectrum  is  made  using  estimates  for  the  CaO  concentration  and 
temperature.  The  CaO  radiator  is  expected  to  be  in  the  gas  phase  boundary  layer  of  the  re¬ 
entry  object.  To  first  order,  the  radiator  thickness  is  not  important.  The  temperature  is  used  to 
compute  the  density  of  the  radiating  state  from  the  CaO  concentration.  The  radiation  is  then 
computed  from  the  usual  equation  for  molecular  radiation  from  diatomic  gases.*’  ’ 

The  computation  of  the  spectrum  requires  the  Einstein  coefficient  (inverse  lifetime)  of 
the  each  of  the  vibrational  transitions.  These  lifetimes  have  been  measured  for  many 
molecules  but  not  for  CaO.  There  are  only  limited  measurements  and  calculation  of  some  of 
the  transitions.  The  Einstein  coefficients  will  then  be  estimated  by  computation  of  Franck- 
Condon  factors  in  the  usual  manner.'"  The  Einstein  coefficient  is  proportional  to  the  cube  of 
the  transition  frequency  times  the  Franck-Condon  factor  times  the  electronic  transition 
moment.  The  electronic  transition  moment  must  be  estimated. 
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Figure  1 .  Emission  Spectrum  of  Interest 


The  computation  of  the  Franck-Condon  factors  requires  information  on  the  potential 
functions.  These  are  found  from  the  spectroscopic  constants  for  the  molecular  states.  The 
spectroscopic  constants  used  are  given  in  Table  1 .  The  potential  functions  are  then  used  in  the 
Schrodinger  equation  to  find  the  vibrational  transition  moments  and  thus  the  Franck-Condon 
factors  and  r-centroids. 


Table  1.  Spectroscopic  Constants'’  for  CaO 


Te 

we 

wexe 

weye 

Be 

alpha 

De 

DO 

A 

11554.8 

718.9 

2.11 

0.40592 

0.00137 

5.4E-7 

37400 

X 

0 

732.03 

4.8307 

0.00887 

0.44444 

0.003282 

6.5E-7 

48350 

The  Franck-Condon  factors  and  r-centroids  were  computed  with  the  assistance  of  the 
codes  provided  by  Le  Roy.  The  Franck-Condon  factors  and  r-centroids  were  found  to  be 
acceptably  close  to  limited  available  literature  values. 

The  values  of  the  electronic  transition  moments  were  then  evaluated  by  comparison  with 
available  lifetime  estimates.  Irwin  and  Dagdigian'^  estimated  the  lifetime  of  the  A'?v'=6) 
state  as  1 55±60  ns.  Diffenderfer  and  Yarkanoy'*  computed  this  lifetime  to  be  130  ns.  Plane 
and  Nien'^  measured  the  lifetime  of  the  radiation  from  the  A'?v'<2)  states  as  150  ns.  Thus, 
we  have  taken  the  sum  of  the  Einstein  coefficients  of  each  of  the  upper  vibrational  levels  to  be 
equal  to  7  ps"'  (each  of  these  states  having  the  same  electronic  transition  moment).  The 
individual  values  are  given  in  Table  2. 

There  is  a  reference'^  to  an  unpublished  evaluation  of  the  variation  of  this  electronic 
transition  moment  with  r-centroid,  constant*exp(-3.58*rv-,v”)-  However,  this  did  not  give 
results  consistent  with  the  lifetime  measurements. 
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Table  2.  Einstein  Coefficients  for  A'i(v')— >  X’lfv"),  jas'' 


vA  V*’ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

3.504 

2.680 

0.734 

0.080 

0.003 

I 

3.125 

0.141 

2.213 

1.302 

0.210 

0.008 

2 

1.717 

1.820 

0.184 

1.318 

1.593 

0.352 

0.015 

3 

0.722 

2.172 

0.602 

0.654 

0.664 

1.692 

0.476 

0.019 

4 

0.254 

1.425 

1.772 

0.080 

0.903 

0.291 

1.677 

0.578 

0.021 

5 

0.086 

0.680 

1.761 

1.143 

0.009 

0.920 

0.104 

1.625 

0.651 

0.018 

0.003 

6 

0.022 

0.273 

1.II2 

1.739 

0.614 

0.138 

0.800 

0.028 

1.561 

0.698 

0.015 

7 

0.005 

0.099 

0.531 

1.421 

I.5I2 

0.259 

0.299 

0.657 

0.006 

1.526 

0.686 

8 

0.002 

0.030 

0.257 

0.967 

1.854 

1.404 

0.084 

0.049 

0.590 

1.762 

9 

0.009 

0.099 

0.541 

1.541 

2.287 

1.257 

0.010 

0.733 

0.522 

10 

0.002 

0.026 

0.196 

0.792 

1.816 

2.133 

0.876 

0.007 

0.783 

0.368 

90  9 1 

The  spectral  radiation  was  computed  ’  for  a  re-entry  object  at  25  km,  observed  at  10 
km  from  a  distance  of  about  250  km.  The  surface  temperature  was  taken  as  1750  K  and  the 
concentration  and  temperature  of  CaO  varied.  The  result  is  shown  in  Figure  2  for  a 
temperature  of  2500K,  a  Na  concentration  of  1  ppm,  and  a  CaO  concentration  of  20  ppm. 

The  proper  evaluation  of  this  radiation  requires  detailed  object  definition,  trajectory, 
ablation  and  flow-field  calculations.  The  purpose  here  is  to  show  document  the  procedure,  to 
the  result  and  to  demonstrate  the  capability. 


wav*  length  (angstroms) 


Figure  2.  Spectral  Radiation  for  a  Reentry  Object  at  25km 


3.  SPACE  CHEMISTRY 

The  application  of  chemical  techniques  to  the  development  of  new  technologies  for  Air 
Force  space-based  systems  and  related  issues  were  also  addressed  under  this  contract.  These 
techniques  were  pursued  in  the  prime  interest  of  enabling  space-based  systems  to  be 
miniaturized  down  to  the  pico-satellite  regime  (<lkg). 

Research  efforts  in  support  of  this  topic  were  towards  the  development  of  a 
miniaturized,  satellite  skin  sensor  technology  capable  of  measuring  ionospheric  neutral 
particle  densities.  This  real-time,  in-situ  capability  enables  the  concept  of  distributed,  remote 
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density  measurements,  via  pico-satellites,  to  provide  3-dimensional  mapping  of  ionospheric 
phenomena  relevant  to  space  weather  forecasting  software  such  as  AF-GEOSpace.  In 
addition,  for  distributed  pico-satellite  systems  that  require  precision  formation  flying,  this 
same  sensor  technology  would  enable  on-line  orbital  drag  evaluation,  to  allow  for  accurate 
station  keeping. 

All  aspects  of  the  research  were  supported  including  data  acquisition,  data  analysis, 
equipment  maintenance  and  upgrade  designs.  All  data  sets  were  delivered  to  AFRL  as  they 
became  available.  In  addition,  day-to-day  operations  of  the  AFRL/RVBXT  space  chemistry 
laboratory  were  addressed  under  this  contract.  This  included  the  oversight  of  the  chemical 
inventory  and  hazardous  waste  accumulation.  Results  of  our  research  were  presented  at 
national  and  international  conferences  and  published  in  peer  reviewed  journals.  Section  6 
lists  five  publications  that  resulted  from  our  efforts  under  this  contract. 

Some  of  the  more  notable  research  results  are  summarized  in  the  following  sections. 

The  investigations  described  were  conducted  at  the  AFRL/RVBXT  facilities.  These  facilities 
include  a  state-of-the-art  Space  Chemistry  Lab  guided-ion  beam  instrument,  and  a  similar 
modular  system  designed  for  simple  modifications  of  the  basic  ion  beam  experiment. 

3.1.  A  Guided-Ion  Beam  Study  of  the  O^S)  +  NH3  System 

We  measured  absolute  cross  section  for  the  reaction  of  ground-state  with  ammonia  at 
collision  energies  in  the  range  from  near-thermal  to  approximately  1 5  eV,  using  the  guided- 
ion  beam  (GIB)  method.  Measurements  were  also  performed  using  ammonia-i/s  to  aid  in  mass 
assignments.  The  reaction  is  dominated  at  low  collision  energies  by  charge  transfer;  however, 
the  cross  section  for  this  exothermic  channel  is  rather  small,  decreasing  sharply  with  energy 
Ifom  ~40  for  normal  ammonia  at  near-thermal  energies  and  leveling  off  at  3.7  above  6 
eV;  the  cross  section  is  slightly  smaller  for  ammonia-i/s.  Other  channels,  corresponding  to  the 
production  of  NH2^  and  NO^,  and  possibly  OH^,  were  detected.  The  NO^  channel,  although 
nominally  exothermic,  is  very  small  and  exhibits  a  threshold  at  ~7  eV.  Product  recoil  velocity 
distributions  were  also  determined  at  selected  collision  energies,  using  GIB  time-of-flight 
methods. 

3.2.  Reactions  of  with  CnH2n+25  n  =  2-4:  A  Guided-ion  Beam  Study 

We  measured  absolute  reaction  cross  sections  for  the  interaction  of  01  with  ethane, 
propane,  and  «-butane  at  collision  energies  in  the  range  from  near  thermal  to  approximately 
20  eV,  using  the  guided-ion  beam  (GIB)  technique.  We  have  also  measured  product  recoil 
velocity  distributions  using  the  GIB  time-of-flight  (TOF)  technique  for  several  product  ions  at 
a  series  of  collision  energies.  The  total  cross  sections  for  each  alkane  are  in  excess  of  100  A2 
at  energies  below  ~2  eV,  and  in  each  case  several  ionic  products  arise.  The  large  cross 
sections  suggest  reactions  that  are  dominated  by  large  impact  parameter  collisions,  as  is 
consistent  with  a  scenario  in  which  the  many  products  derive  from  a  near-resonant, 
dissociative  charge-transfer  process  that  leads  to  several  fragmentation  pathways.  The  recoil 
velocities,  which  indicate  product  ions  with  largely  thermal  velocity  distributions,  support  this 
picture.  Several  product  ions,  most  notably  the  C2H3-I  fragment  for  each  of  the  alkanes, 
exhibit  enhanced  reaction  efficiency  as  collision  energy  increases,  which  can  be  largely 
attributed  to  endothermic  channels  within  the  dissociative  charge  transfer  mechanism. 
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3.3.  Hyperthermal  Reactions  of  0^(^8372)  with  CD4  and  CH4 


Theoretical  and  experimental  methods  were  applied  to  the  study  of  the  reaction 
dynamics  in  hyperthermal  collisions  of  (^83/2)  with  methane.  Measurements  of  the  absolute 
reaction  cross  sections  for  the  interaction  of  with  CD4  and  CH4  were  obtained  at  collision 
energies  in  the  range  from  near-thermal  to  approximately  15  eV,  using  the  guided-ion  beam 
(GIB)  technique.  Product  recoil  velocity  distributions,  using  the  GIB  time-of-flight  (TOP) 
methods,  were  determined  for  several  product  ions  at  selected  collision  energies.  The  main 
reaction  channel,  charge  transfer,  proceeds  via  large  impact  parameter  collisions.  A  number  of 
minor  channels,  involving  more  intimate  collisions,  were  also  detected.  Ab  initio  electronic 
structure  calculations  have  been  performed  with  different  levels  of  theory  and  basis  sets, 
including  high-level  coupled-cluster  calculations  to  determine  the  energies  of  reaction 
intermediates  and  transition  states  for  reaction.  Several  reaction  paths  on  both  quartet  and 
doublet  electronic  states  of  (O  CH4)^  are  found,  and  these  provide  a  reasonable  qualitative 
interpretation  of  the  experiments.  Although  most  of  the  products  can  be  produced  via  spin- 
allowed  pathways,  the  appearance  of  CHs^  at  low  energies  suggests  that  intersystem  crossing 
plays  some  role. 


4.  SMEI  SUPPORT  AND  INTERPLANETARY  RESEARCH 

One  of  the  goals  of  our  research  efforts  on  this  contract  was  to  work  toward 
understanding  the  origins  and  interplanetary  propagation  (IP)  of  disturbances  that  affect  the 
geosphere,  and  to  support  and  analyze  the  data  from  the  AFRL/NASA  Solar  Mass  Ejection 
Imager  (SMEI)  experiment  which  was  launched  on  the  Coriolis  spacecraft  in  January  2003. 
Our  research  utilized  pertinent  satellite  and  ground-observations  of  solar,  IP  and  geomagnetic 
phenomena  and,  in  particular,  the  new  observations  from  the  unique  SMEI  experiment. 

Our  studies  involved  both  historical  and  new  data  acquired  over  the  course  of  this 
contract.  It  primarily  addressed  three  general  topics  including: 

■  the  solar  eruptive  phenomena  which  lead  to  sporadic  geomagnetic  storms 

■  IP  shocks  and  other  signatures  of  solar  ejecta 

■  the  characteristics  of  IP  disturbances  which  produce  geomagnetic  storms 

The  primary  results  of  our  research  were  presented  in  numerous  conferences  nationally 
and  internationally.  Two  papers  directly  related  to  the  work  performed  under  this  contract 
were  published  in  peer-reviewed  journals.  They  are  included  in  Section  6.0.  The  most 
relevant  research  efforts  are  summarized  in  the  following  sections. 

4.1.  Major  Geomagnetic  Storms 

Seventy- nine  major  geomagnetic  storms  (minimum  Dst  <-100  nT)  observed  in  1 996 
to  2004  were  the  focus  of  a  “Living  with  a  Star”  Coordinated  Data  Analysis  Workshop 
(CDAW)  in  March  2005.  In  nine  cases,  the  storm  driver  appears  to  have  been  purely  a 
corotating  interaction  region  (CIR)  without  any  contribution  from  coronal  mass  ejection- 
related  material  (interplanetary  coronal  mass  ejections  (ICMEs)).  These  storms  were 
generated  by  structures  within  CIRs  located  both  before  and/or  after  the  stream  interface  that 
included  persistently  southward  magnetic  fields  for  intervals  of  several  hours.  We  compare 
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their  geomagnetic  effects  with  those  of  159  CIRs  observed  during  1996-2005.  The  major 
storms  form  the  extreme  tail  of  a  continuous  distribution  of  CIR  geoeffectiveness  which  peaks 
at  Dst  =  -40  nT  but  is  subject  to  a  prominent  seasonal  variation  of  ±40  nT  which  is  ordered 
by  the  spring  and  fall  equinoxes  and  the  solar  wind  magnetic  field  direction  toward  or  away 
from  the  Sun.  The  O'Brien  and  McPherron  (2000)  equations,  which  estimate  Dst  by 
integrating  the  incident  solar  wind  electric  field  and  incorporating  a  ring  current  loss  term, 
largely  account  for  the  variation  in  storm  size.  They  tend  to  underestimate  the  size  of  the 
larger  CIR-associated  storms  by  Dst  =  20  nT.  This  suggests  that  injection  into  the  ring  current 
may  be  more  efficient  than  expected  in  such  storms.  Four  of  the  nine  major  storms  in  1996- 
2004  occurred  during  a  period  of  less  than  three  solar  rotations  in  September  to  November 
2002,  also  the  time  of  maximum  mean  IMF  and  solar  magnetic  field  intensity  during  the 
current  solar  cycle.  The  maximum  ClR-storm  strength  found  in  our  sample  of  events,  plus 
additional  23  probable  CIR-associated  Dst  <  -100  nT  storms  in  1972-1995,  is  {Dst  =  -161 
nT).  This  is  consistent  with  the  maximum  storm  strength  {Dst  =  -180  nT)  expected  from  the 
O'Brien  and  McPherron  equations  for  the  typical  range  of  solar  wind  electric  fields  associated 
with  CIRs.  This  suggests  that  CIRs  alone  are  unlikely  to  generate  geomagnetic  storms  that 
exceed  these  levels. 

4.2.  Solar  Mass  Ejection  Imager  (SMEI)  Observations  of  CMEs 

The  Solar  Mass  Ejection  Imager  (SMEI)  on  the  Coriolis  spacecraft  has  been  obtaining 
white  light  images  of  nearly  the  full  sky  every  102  minutes  for  three  years.  We  present 
statistical  results  of  analysis  of  the  SMEI  observations  of  coronal  mass  ejections  (CMEs) 
traveling  through  the  inner  heliosphere;  139  CMEs  were  observed  during  the  first  1.5  years  of 
operations.  At  least  30  of  these  CMEs  were  observed  by  SMEI  to  propagate  out  to  1  AU  and 
beyond  and  were  associated  with  major  geomagnetic  storms  at  Earth.  Most  of  these  were 
observed  as  frontside  halo  events  by  the  SOHO  LASCO  coronagraphs. 


5.  DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM  (DMSP) 

Our  efforts  in  support  of  DMSP  research  included  the  development  of  a  time  history 
data  base.  This  data  base  includes  measurements  of  the  early  DMSP  vehicles  F13,  F14  and 
FI  5  together  with  the  measurements  made  by  DMSP  vehicles  FI 6  and  FI 7,  launched  during 
the  course  of  this  contract.  These  databases  have  been  essential  to  the  generation  of  data 
products  of  AFRL/RVBX  use  as  well  as  space  weather  applications  at  other  government 
facilities  such  as  the  Air  Force  Weather  Agency  (AFWA). 

5.1.  DMSP  Data  Base  Development 

The  objective  of  this  effort  was  to  maintain  the  time  history  databases  from  the 
previously  flown  DMSP  vehicles  and  to  continue  to  develop  time  history  databases  for  the 
vehicles  launched  during  the  course  of  this  contract. 

To  date,  the  data  base  includes  measurements  from  the  DMSP  vehicles,  F13,  F15,  F16 
and  FI 7.  AFWA  reported  a  hardware  failure  on  the  F14  recorded  on  August  23,  2008.  Thus 
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no  data  was  received  for  FI  4  after  this  date.  The  tasks  inherent  in  the  development  of  this 
database  are  described  in  this  section  of  the  report. 

AFWA  redesigned  the  DMSP  mission  sensor  satellite  data  files,  known  as  PREP  files, 
into  Raw  Sensor  Data  Record  (RSDR)  files.  The  new  format  replaced  the  36-bit  processed 
files,  containing  sensor  and  file  structure  specific  information,  with  a  standard  32-bit  format 
containing  only  sensor  specific  information.  Boston  College  efforts  included  modifying  pre¬ 
developed  software  to  transform  the  RSDR  files  into  time-ordered  databases.  This  software 
reads  the  RSDR  files  and  generates  SSJ/4,  SSIES  and  SSM  time-ordered  databases.  Boston 
College  personnel  run  the  software  and  provide  quality  control  checks  on  a  daily  basis. 

Each  RSDR  file  contains  one  payload  of  sensor  data.  Instead  of  receiving  twelve  hours 
of  data  bundled,  each  payload  is  received  soon  after  it  is  processed.  Calculated  ephemeris 
data,  time  codes,  and  on-board  ephemeris  are  provided  for  each  second  of  data  rather  than  for 
each  minute.  This  time  and  ephemeris  information  is  identical  for  specific  DMSP  flights. 

Each  36-bit  sensor  data  word  is  divided  into  three,  12-bit  words  and  stored  right  justified  in  a 
16-bit  word.  However,  the  RSDR  file  does  not  contain  the  year  and  day  associated  with  the 
sensor  data.  The  year  and  day  must  be  calculated  by  the  use  of  a  set  of  subroutines  that 
assume  that  a  given  time  and  its  associated  ephemeris  are  accurate  and  that  the  day  given  is 
the  day  of  data  or  a  day  within  the  past  three  or  four  days.  Each  RSDR  file  will  have  a  name 
to  uniquely  identify  the  data  set.  The  file  name  will  contain  the  satellite  designation,  playback 
revolution  number,  the  year,  day,  hour,  and  minute  that  the  RSDR  was  created,  and  the  sensor 
identification. 

The  new  Phase  1  time  history  database  formats  contain  only  16-bit  or  32-bit  integer  data 
and  8-bit  ASCII  characters.  Each  9-bit  data  sensor  value  is  stored  right  justified  in  a  1 6-bit 
word.  Each  block  of  data  contains  one  minute  of  data. 

DMSP  data,  stored  in  RSDR  format,  is  transformed  into  time-ordered  databases  on  a 
UNIX  workstation.  Raw  data  sets  for  the  SSJ,  SSIES,  and  SSM  are  acquired,  copied,  and 
archived. 

Input  data  is  copied  by  using  Secure  File  Transfer  Protocol  (SFTP)  between  the 
computer  at  Offutt  AFB  and  the  Amber,  a  SUN  Microsystems  UNIX  computer.  To  automate 
this  procedure  a  script  is  executed  at  designated  times  during  the  day  to  acquire  these  files. 

For  each  experiment,  for  each  vehicle,  a  multi-day  period  is  stored  on  UniTree  on  a  daily 
basis  and  is  also  backed  up  on  DVD’s  on  a  monthly  basis. 

Computer  modules,  written  in  FORTRAN,  with  the  input  data  in  RSDR  format,  were 
carefully  developed  to  generate  SSJ,  SSIES,  and  SSM  time-ordered  databases.  Data  values 
are  stored  in  raw  counts  for  each  second;  ephemeris  values  are  stored  at  each  even  minute.  A 
common  software  package  that  computes  geomagnetic  parameters  is  used  for  each 
experiment.  A  subroutine  was  written  to  unpack  the  three,  12-bit  “data”  words  and  convert 
them  to  four,  9-bit  sensor  words.  Since  any  one  revolution  will  neither  start  nor  end  exactly  at 
the  start  or  end  of  a  day,  care  must  be  taken  to  ensure  that  the  set  of  files  used  as  input  for  one 
day  overlap  into  the  preceding  and  succeeding  days. 

To  facilitate  the  creation  of  a  database,  script  files  were  written  for  each  experiment  that 
generate  a  procedure  that  process  multiple  days  of  data  that  consist  of  multiple  input  files.  As 
the  script  procedure  is  executed  a  log  file  is  created  that  contains  processing  times  and  error 
indications.  Directories  for  each  experiment  are  created.  A  listing  of  any  directory  contains 
the  name,  size,  and  date  created  of  each  individual  file  processed.  If  the  file  size  for  a  day  of 
data  is  not  what  is  expected  for  a  particular  experiment,  it  will  be  verified  that  the  input  data  is 
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unavailable  or  the  error  will  be  found  and  the  day  processed  again.  The  database  for  each  day 
is  stored  on  the  UniTree.  Copies  of  the  SSJ  database  files  are  sent  on  a  daily  basis  via  FTP  to 
Three  outside  agencies.  Copies  of  the  SSIES,  SSJ  and  SSM  database  files  are  also  archived  on 
to  DVD’s  on  a  monthly  basis. 

5. 1. 1.  RSDR  SSIES 2  Phase  II  Processing 

The  Phase  II  processing  system  accesses  the  Phase  I  time  history  database.  This  system 
is  comprised  of  several  software  routines  coded  in  FORTRAN  and  is  driven  by  a  UNIX 
script.  Telemetry  data  values  from  the  Driftmeter  (DM),  Scintillation  Monitor  (SM), 

Retarding  Potential  Analyzer  (RPA),  Langmuir  Probe  (LP)  and  Microprocessor  (MP) 
experiments  are  converted  into  geophysical  parameters  and  written  to  output  files  at  one- 
minute  intervals.  The  automated  processing  is  done  on  a  daily  basis  with  a  crontab  job 
initiating  the  shell  script  at  a  specified  time  in  the  morning.  The  files  are  archived  on  the 
UniTree  file  system.  Copies  of  the  Phase  II  data  base  files  are  sent  on  a  daily  basis  via  FTP  to 
two  outside  agencies. 

5. 1.2.  RSDR  SSIES2  and  SSIES3  Phase  II EDR  Processing 

The  processing  consists  of  accessing  the  phase  I  files,  converting  the  files  to  RSDR 
format,  then  to  EDR  data  files.  The  EDR  files  are  compressed  and  then  archived  to  the 
UniTree  file  system. 

A  copy  of  the  compressed  file  is  also  transferred  to  the  FTP  server  (dropbox)  so  that  the 
Naval  Research  Laboratory  can  retrieve  the  files  via  FTP  and  to  FTP  servers  at  NGDC  and 
the  Aerospace  Corporation. 

The  processing  is  done  by  a  combination  of  UNIX  shell  scripts,  Perl  scripts  and 
FORTRAN  programs.  The  automated  processing  is  done  on  a  daily  basis  with  a  crontab  job 
initiating  the  shell  script  at  a  specified  time  in  the  morning. 

5. 1. 3.  RSDR  SSM  Phase  II MFR  Processing 

The  processing  consists  of  accessing  the  phase  I  files,  converting  the  files  to  RSDR 
format,  then  to  MFR  (ASCII)  data  files.  Another  program  is  run  converting  the  ASCII  files  to 
binary,  the  binary  files  are  compressed  then  transferred  to  FTP  servers  at  JHU/APL,  NGDC 
and  the  Aerospace  Corporation.  Both  types  of  files  are  also  archived  to  the  Unitree  file 
system.  The  processing  is  done  by  a  combination  of  UNIX  shell  scripts,  Perl  scripts  and 
FORTRAN  programs.  The  automated  processing  is  done  on  a  daily  basis  with  a  crontab  job 
initiating  the  shell  script  at  a  specified  time  in  the  morning. 

5. 1. 4.  DMSP3  and  SWX  Web  Survey  Plot  Files  (SSIES,  SSJ  &  SSM) 

The  DMSP3  processing  consists  of  a  shell  script  that  accesses  input  files  from  Unitree, 
then  runs  IDL  programs  that  generate  High  and  Low  Latitude  data  survey  plot  files.  The  plots 
files  are  copied  to  the  Web  server  for  users  to  access  inside  the  firewall.  The  automated 
processing  is  done  on  a  daily  basis  with  a  crontab  job  initiating  the  shell  script  at  a  specified 
time  in  the  morning.  The  High  latitude  web  plots  include  SSJ  Ion  and  Electron;  SSIES  Ni  and 
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N  [0+],  Horizontal  and  vertical  drift  of  thermal  ions,  and  SSM  measured  minus  model 
magnetic  field  data.  The  Low  latitude  web  plots  include  SSIES  Ni/N[0+],  Ion  Drift  meter, 
Temperature  and  SM  total  Ion  density  survey  data.  Another  script,  also  initiated  by  a  crontab 
job  is  run  on  SWX  that  copies  the  plots  files  generated  on  dmsp3  on  to  the  SWX  web  server 
with  a  two  day  delay  for  users  to  access  outside  of  the  firewall.  The  log  files  and  web  pages 
are  viewed  on  a  daily  basis  to  check  for  errors  in  processing 

5. 1. 5.  Raw  Sensor  Data  Record  Files  (RSDR) 

RSDR  files  for  SSM,  SSIES  and  SSJ  for  the  period  of  August  1  through  October  3 1 , 
2008  were  acquired  from  AFWA  for  vehicles  F13,  F15,  F16  and  F17  (The  F14  recorder  failed 
on  August  23  ).  AFWA  had  numerous  problems  during  this  reporting  period  including 
missing  files,  corrupted  data,  especially  at  the  transition  of  one  day  to  the  next  and  repeating 
rev  readouts  that  also  contained  bad  data.  The  problems  were  reported  to  AFWA.  Many  days 
included  some  missing  files.  In  some  cases  they  were  able  to  recover  the  missing  files.  When 
missing  files  were  recovered,  all  of  the  follow-on  programs  for  the  affected  data  were  rerun. 

SSJ,  SSIES  and  SSM  Time  History  database  files  which  use  RSDR  files  as  input  have 
been  generated  for  the  period  of  August  1  through  October  31,  2008.  SSIES  Phase  II  database 
files  which  use  THDB  files  as  input  have  also  been  generated  for  the  above  time  period  for 
vehicles  F 13  and  FI 5. 

SSIES2/3  FI 3,  FI  5,  FI 6  and  FI 7  EDR  Phase  II  files  were  processed  from  August  1 
through  October  3 1 ,  2008.  A  request  by  the  program  manager  to  run  the  SSIES2/3  Phase  II 
EDR  program  with  the  option  to  create  24  samples  per  second  files  of  SM  output  was 
completed  for  the  period  of  July  14-28  2008.  A  UNIX  script  was  written  to  process  the  data 
further  listing  only  time  and  density  values. 

SSM  Phase  II  MFR  processing  was  completed  from  August  1  -  October  31  2008. 

SSIES,  SSJ  and  SSM  Time  History  database  files,  SSIES  Phase  II  database  files,  SSIES  EDR, 
SSM  MFR  and  Ephemeris  files  have  been  archived  onto  DVD’s  from  July  1  -  September  30 
2008.  Two  copies  of  the  archive  DVD’s  were  made  and  sent  to  outside  agencies. 

RSDR  raw  files  for  SSJ,  SSIES  and  SSM  were  archived  onto  DVD’s  for  the  period  of 
July  1  through  September  30,  2008. 

DMSP3  Web  server  plot  files  were  generated  for  the  period  of  August  1  through 
October  3 1 .  The  plot  files  were  copied  to  the  SWX  web  server  through  October  29. 

An  anomaly  was  reported  to  the  program  manager  pertaining  to  the  FI  5  Low  Lat  SM  ION 
density  survey  plots.  The  problem  started  on  September  10.  It  was  determined  that  that  the 
instrument  may  need  to  be  reset.  Also  the  programmer  responsible  for  the  SSIES  EDR  codes 
made  updates  to  the  codes  so  the  FI 5  EDR  data  from  September  to  current  will  need  to  be 
rerun. 

5.2.  DMSP  J4/J5  Data  Management 

Boston  College  maintained  a  complete  archive  of  all  DMSP  J  data  from  1983  through 
the  present,  beginning  with  the  F6  data  and  continuous  through  the  currently  active  DMSP 
satellites  (F13,  F15,  F16  and  F17).  This  archive  is  maintained  with  triple  redundancy  on  a 
central  AFRL/RV  file  system  with  a  fourth  copy  of  CD.  This  archive  has  been  made 
available  to  all  AFRL  personnel.  In  addition,  several  levels  of  processed  summary  data  files 


13 


are  maintained  along  with  viewing  software  in  support  of  various  AFRL  projects.  Boston 
College  personnel  use  the  archive,  processed  files  and  software  tools  to  monitor  the  status  of  J 
instruments  on  currently  active  DMSP  satellites. 

The  archive  process  has  been  fully  automated,  but  requires  manual  processing  whenever 
network  problems  interrupt  the  automation  string,  which  occurs  on  most  weekends.  In 
addition,  log  files  were  routinely  monitored  for  quality  control. 

During  this  contract  period,  the  processing  software  has  undergone  several  updates  and 
viewing  software  has  been  upgraded  for  more  rapid  use  on  newer  computers. 

5.3.  DMSP  J5  Calibration 

Boston  College  personnel  supported  the  calibration  and  delivery  the  J5  instruments  for 
the  DMSP  program.  Calibrations  were  performed  using  AFRL’s  calibration  facility  for 
determining  detector  response  to  ions  and  electrons  with  energies  between  1 00  eV  and  20 
KeV.  The  first  J5  calibration  using  this  facility  was  completed  in  June  2000  and  the  second  in 
August  2002. 

Several  serious  problems  and  many  minor  problems  with  the  calibration  facility 
hardware  and  software  were  detected  during  the  initial  J5  calibration.  Since  then,  new  motors 
and  motor  controllers  with  standard  drivers  were  installed  and  a  Labview  user  interface  was 
written  to  replace  the  existing  obsolete  calibration  software. 

Calibration  results  for  the  S20  and  S16  and  SI 7  were  completed  and  released  to  the  Air 
Force  DMSP  program. 

Software  for  the  calibration  system,  CalSys4,  was  continuously  upgraded,  debugged, 
tested  and  documented.  By  the  end  of  this  contract,  a  draft  user’s  guide  was  developed  and 
provided  to  AFRL. 

5.4.  DMSP  Auroral  Boundary  Maintenance 

Boston  College  researchers  have  maintained  a  random  access  database  of  auroral 
boundaries  determined  from  DMSP  J  data  that  is  complete  from  1982  to  the  present.  The 
algorithm  for  determining  the  auroral  boundaries  was  developed  by  Boston  College 
personnel.  During  the  course  of  this  contract,  we  have  continued  to  maintain  and  upgrade  all 
relevant  software.  In  addition,  we  have  been  provided  the  algorithm  and  operating  system  to 
the  Air  Weather  Service.  The  boundary  data  has  been  used  as  an  indicator  for  solar  activity  in 
their  Space  Weather  Forecasting  System. 

The  bounding  algorithm  requires  six  threshold  parameters  that  must  be  determined  for 
each  new  DMSP  satellite  to  produce  results  consistent  with  the  archived  boundary  data.  In 
addition,  yearly  updates  to  the  data  files  have  provided.  At  the  close  of  this  contract,  routine 
data  processing  and  data  quality  control  was  up  to  date.  In  addition,  as  set  of  threshold 
parameters  for  the  AWS  boundary  algorithm  for  F16  and  FI  7  was  delivered  as  soon  as 
possible  following  the  successful  launch  of  these  vehicles. 

5.5.  CEASE  land  11 

The  Compact  Environmental  Anomaly  Sensor  Experiment  (CEASE)  consisted  of  a 
small  environmental  hazard  detection  system  that  was  both  lightweight  and  had  low  power 
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requirements.  The  total  package  provided  data  on  radiation  dose,  single  event  upsets,  and 
charging  (spacecraft  and  dielectric). 

A  software  system  has  been  developed  that  will  take  the  prime  data  set,  the  CEASE 
Science  file,  and  parse  the  data  into  individual  time  tagged  data  base  files.  A  program  to 
create  an  ephemeris  time  history  database  from  a  reconstructed  ephemeris  file  has  been 
written.  CEASE  was  being  flown  on  the  Tri-Service  Experiments  Mission  5  (TSX-5)  vehicle 
that  stopped  generating  data  on  July  5,  2006.  Since  the  mission  was  initiated  in  2000,  the  data 
set  is  comprised  of  six  years  of  continuous  coverage.  Figure  3  illustrates  the  survey  of  the  >  1 
MeV  Electrons  in  an  L-shell  versus  Day  for  the  whole  missions. 


VEAR 

Figure  3.  CEASE  Measurements  from  the  TSX-5  Vehicle  (July  2000-July  2006) 

The  objective  of  this  project  was  to  generate  and  maintain  the  CEASE  science  and 
ephemeris  databases.  This  has  been  achieved  and  the  following  describes  the  contents  of  the 
database: 

■  For  the  TSX-5  vehicle:  CEASE  science  and  ephemeris  databases  have  been  generated 
and  cover  the  time  period  of  Day  160,  Year  2000  to  Day  186,  Year  2006 

■  For  the  DSP-21  vehicle,  second  vehicle  equipped  with  a  CEASE  sensor:  CEASE 
science  and  ephemeris  databases  have  been  generated  and  cover  the  time  period  of  Day 
221,  Year  2001  to  Day  117,  2009. 
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